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ABSTRACT 

We present a photometric catalog of 8,735,004 proper motion selected low-mass stars (KML-spectral 
types) within the Sloan Digital Sky Survey (SDSS) footprint, from the combined SDSS Data Release 10 
(DRIO), Two-Micron All-Sky Survey (2MASS) Point Source Catalog (PSC), and Wide-field Infrared 
Survey Explorer {WISE) AllWISE catalog. Stars were selected using r — i, i — z, r — z, z — J, and 
z — W1 colors, and SDSS, WISE, and 2MASS astrometry was combined to compute proper motions. 
The resulting 3,518,150 stars were augmented with proper motions for 5,216,854 earlier type stars 
from the combined SDSS and United States Naval Observatory Bl.O catalog (USNO-B). We used 
SDSS-I-USNO-B proper motions to determine the best criteria for selecting a clean sample of stars. 
Only stars whose proper motions were greater than their 2a uncertainty were included. Our Motion 
Verified Red Stars (MoVeRS) catalog is available through SDSS CasJobs and VizieR. 

Keywords: catalogs — infrared: stars — proper motions — stars: low-mass — stars: kinematics and 
dynamics — stars: late-type 


1. INTRODUCTION 

Over the past century, photometric surveys—both dig¬ 
ital and photographic—have played an important role 
in many facets of astronomy. One of the largest limita¬ 
tions to these surveys has been object classification for 
point sources that have similar colors and morphologies 
(e.g., M giants, M dwarfs, QSOs, and distant luminous 
red galaxies). One method for separating nearby stellar 
populations from more distant objects is measuring tan¬ 
gential motion on the sky, or proper motion. Proper mo¬ 
tions are also important for distinguishing and investigat¬ 
ing kinematically distinct populations within our Galaxy 
(e.g., moving groups, disk and halo stars, etc.). 

A number of large, all-sky catalogs of stellar positions 
and proper motions now exist. The United States Naval 
Observatory (USNO) has had a long history of track¬ 
ing astrophysical obje cts, starting with th eir first pub¬ 
lished catalog UJl.O (Monet et ah 19941, an d subse- 

USNO- 


quent ly repli ed by USi\U-Al.(J ( [Monet 


9961 


A2.0 (|^^|^et| |1998 ), and ultimately USiNU-Bl.O (Monet 


et al. |2003 nereatter USNO-B). These catalogs are the 
result ot the Precision Measuring Machine at the USNO 
Flagstaff Station, undertaking a photometric survey over 
^50 years on Schmidt plates. USNO-B culminated in 
a catalog containing the positions, proper motions, and 
magnitudes for over a billion objects. However, the 
proper motions and positions in USNO-B are relative, 
not absolute, making it difficult to compare observations 
with later epoch observation s on a well-defined system 
(e.g.,|R6ser et al.|2008l|2010D. Other proper motion cat¬ 
alogs ombteTncIuHel but are not limited to): the SU- 
PERBLI NK catalog of norther n stars with large proper 
motions (Lepine & Shara 2005 hereafter LSPM), the Po¬ 
sitions and Proper Motions catalog (PPM, PPMX, PP- 

Electronic address: ctheisse@bu.edu 

^ Department of Astronomy, Boston University, 725 Common¬ 
wealth Avenue, Boston, MA 02215, USA 

^ Department of Physical Sciences, Embry-Riddle Aeronautical 
University, 600 South Clyde Morris Blvd., Daytona Beach, FL 
32114, USA 


MXL; Roser & Bastian 1993 Roser et al. 2008 2010), 
and catalog s calibrate^ with the bloan Digi tal bky Pur¬ 
vey (SPSS; [York et al.||20dot an d USNO-B (|Munn et al. 
2004| GouIq & Koilmeier||2U04 ) surveys. 


Une ot the newest surveys from the USNO is the CCD 
Astrograph Catalog (UCAC; Zacharias et al.|[^00|), now 
in its 4*'' (and final) data release (UCAC4; Zacharias 


et al. 20131. While UCAC4 does not contain as many 
objects as USNO-Bl.O (~10® objects), and has a shorter 
time baseline (~6 yrs), it is approx imately five times 


more precise (Zacharias et al. 2000) due to the use of 


CCDs instead of photographic plates. UCAC is also in 
an absolute reference frame, the International Coordi¬ 
nate Reference System (ICRS). The newest undertak¬ 
ing (started in April 2012) by the USNO is the USNO 
Robotic Astrometric Telescope ( URAT), currently in it s 
initial data release with (URATl; Zacharias et al.|2015[ ). 
This catalog is expected to achieve precision astromet¬ 
ric measurements (^10 mas) for 500 million sources, but 
will be r elatively shallow in compar ison to other surveys 
{R « 18; |Zacharias & Gaume||2011|), and it will be many 
years before the deepest all-sky data release (northern 
and southern hemisphere) is available. 

The best example to date of high-precision space based 
astromet ry is the s urvey performed by the Hipparcos 
satellite (ESA 1997), making precise astrometric mea- 
surements (< 10 mas) for millions of stars. The cur¬ 
rent reahzationof Hipparcos data is the Tycho-2 cata¬ 
log (H0g et al. ]2000| containing astrometric information 
for approximately 2.5 million stars. Prior to Tycho-2, 
the Hipparcos catalog ( [Perryman et al.|1997 ), containing 
extremely high precision astrometric measurements for 
118,218 stars, defined the ICRS at optical wavelengths. 
The ICRS is the standard reference frame that modern 
surveys (e.g., SDSS, a nd the future Larg e Synoptic Sur¬ 
vey Telescope or LSST; Ivezic et al. 2008) use to calibrate 
their astrometry, since the majority ot these surveys are 
tied to either Tycho-2 or UCAC. 

Although many proper motion catalogs exist, they are 
typically tied to surveys that are biased towards the blue 
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end of the spectrum (e.g., USNO-B). This makes most 
current proper motion catalogs severely incomplete at 


smaller catalogs (e.g.. 

Lepine & Shara 

2005| 

Peacon & 

Hambly 

2007 

I'aherty et al. 

2DD9 

nt 

is only in recent 


years that an all-sky infrared point-source catalog with 
enough astrometric precision to build a more complete 
proper motion catalog for the reddest point sources has 
become available. Combining observations from all-sky 
and large area surveys taken over the past two decades 
will allow us to compute reliable proper motions for the 
lowest-mass stars. 

The past two decades have seen the emergence of three 
of the most important astronomical surveys for stud¬ 
ies of lo w-mass stars: the Tw o Micron All-Sky Survey 
(2MASS; Skrutskie et al. 2006), SP SS, and the Wide-fiel d 
Infrared Purvey Explorer {WISE\ [Wright et al.||2010|). 
2MASS conducted observations between 1997 and 2001 
in three near-infrared (NIR) bands (J: 1.25 /rm, H: 1.65 
/rm, and Kg-. 2.17 /im). The 2MASS point source cata¬ 
log (PSC) contains over 470 million objects. SPSS con¬ 
ducted visible wavelength observations starting in 2000, 
with some observations as recent as the last five years. 
SPSS observed the sky in five visible wavelength bands 
(ugriz), and Pata Release 10 (PRIO; Ahn et al. 20141 
contains approximately 260 million point sources. WISE 
began observing the entire sky in 4 mid-infrared bands 
(3.4, 4.6, 12, and 22 fim) starting in 2010. A second 
post-cryogenic mission (NEOWISE; Mainzer et al. ]MI] ) 
was carried out at the end of 2010 using the two short- 
est bands, and surveying the entire sky over the course 
of a year. The AllWISE catalog combines both WISE 
missions to create a catalog with enhanced photomet¬ 
ric sensitivity and accuracy, and improved astrometric 
precision above each individual mission’s data products. 
The AllWISE catalog contains over 747 million objects. 
Combining these three surveys provides a time baseline 
of ~10 years. 

The ubiquity of M dwarfs throug hout the Galaxy 


■"70% of the total stellar population; Bochanski et al. 


20101, coupled with the fact that M dwarfs have main- 


sequence lifetimes longer than the curre nt age of the Uni¬ 


verse ("^10^^ yrs; Laughlin et al. 19971, make M dwarfs 


important laboratories tor studying numerous aspects 
of astronomy (e.g.. Galactic and stellar evolution, kine¬ 
matics, etc.). Recent results have also suggested that 
M dwarfs have a strong penchant for building terres - 
trial planets (e.g.. Pressing fc Charbonneau|2013 2015). 
This affinity for creating terrestrial planets, coupled with 
the relative ease for finding terrestrial planets around 
M dwarfs (due to their size ratios and small orbital 
distances), make M dwarfs important hosts for study¬ 
ing Earth-sized planets and habitability throughout the 
Galaxy. 

Many large catalogs of M dwarfs (dMs) currently ex¬ 
ist, incl uding the Palomy/Michigan State Un iversity 
(PMSU; [Reid et al.||1995l ||Hawley et al.||1996|) survey 
("^2400 s pectroscopic dM^ and the SPbb Pata Release 


7 (PR7; Abazajian et al. 


2009|) spectroscopic catalog 


(70,841 spectroscopic dlVls-jWest et al. 2011|). However, 
such catalogs make up only a small traction of the mil¬ 
lions of photometric dMs contain ed within SPSS. For 
example, ||Bochanski et al.j (|2010|) used SPSS PR7 to 
retrieve ~15 million photometrically selected, but not 


proper motion verified, dMs to investigate the mass and 
luminosity functions of the Galactic disk. While most of 
the red point sources in SPSS are M dwarfs rather than 
giants or r ed galaxies for the color and magnitude range 
chosen by Bochanski et al. (2010), proper motions can 
help to s elect bona-fide low-mass stars. 

Gaia ([Perryman et al. 2001) is currently conduct¬ 
ing the largest astrometric survey to date. Gaia is a 
magnitude limit ed survey, the limi ts of which are shal¬ 
lower (r < 20; jlvezic et al. 2012) than the combined 
WTS'E+SPSS limits (r < 22.2), making Gaia incomplete 
for the faintest and lowest-mass (reddest) stars. Approx¬ 
imately 60% of the stars in the combined photometric 
dataset of WISE, SPSS, and 2MASS have r > 20, mak¬ 
ing the majority of low-mass stellar candidates beyond 
the reach of Gaia (for relative point-s ource densities at 
different r magnitudes see Figure 6 of [Bochanski et al.j 

Mo|. 

To make use of a larger photometric sample of dMs 
within SPSS, we combine WISE, SPSS, and 2MASS 
observations to compute proper motions over ~10 year 
baselines for photometrically selected objects with dM 
colors. After a color cut and selection of reliable proper 
motions, we are left with a sample of 8,735,004 stars. In 
Section [^ we outline our methods for computing proper 
motions and errors. We also estimate the intrinsic uncer¬ 
tainty within our catalog using SPSS selected quasars. In 
Section [^ we discuss the selection criteria used to build 
our phot ome tric dM sample, and address contamination 
(Section 3.4[ ). In, Section we assess the reliability of 
our catalog. We augment our derived proper motions 
with measurements from SPSS-I-USNO-B in Section [^ 
In Section we discuss the properties of our Motion 
Verified Red Stars (MoVeRS) catalog and how to query 
it. We describe preliminary science results that can be 
achieved with this catalog in Section [^ Our summary 
follows in Section [8| 


2. METHODS: COMBINING WISE, SDSS, AND 2MASS 
2.1. Astrometric Algorithms 

SPSS was originally ca librated against UGAG and 
Tycho-2 |Pier et al.[[2003[), but as of its seventh data 
release (PR7) was calibrated against UCAC2 and an in¬ 
ternal UCAC release known as “rl4^.” This recalibration 
reduced systematic errors from ~75 mas (Tycho-2) and 
"^45 mas (UCAC) to less than 20 mas^. SPSS is also 
on the ICRS since the Tycho-2 catalog is based upon 
Hipparcos astrometry, which defines the ICRS at visible 
wavelengths. All SPSS astrometric calibrations for this 
study were performed in the r-band. 

We noticed that SPSS positional errors (raErr, de- 
cErr) in all SPSS Pata Releases more recent than PR7 
are pixel-centroiding errors rather than absolute astro- 
metric errors (B. A. Weaver, personal communication). 
To compute absolute astrometric errors, we found the 
total error (centroiding plus calibration) in great circle 
coordinates using the following equations. 


(rowcErr X 0.3961)2 _|_ muErr^ (1) 

https://www.sdss3.orK/drlO/algorithms/astrometry.php 
http: //classic. sdss. org/dr7/algorithins/astroinetry. html 
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and 

cr^ = (colcErr X 0.3961)2 _|_ nuErr^, (2) 

where rowcErr and COLCErr are the row center 
and column enter position errors in r—band coordi¬ 
nates, respectively, and are found in the SDSS CasJobs 
“PhotoObj” table. The fields muErr and NuErr are 
the astrometric errors in Great Circle coordinates (/i and 
v) for the r-band, and are found in the “Field” table. The 
factor of 0.3961 is the SDSS pixel scale (arcsec pix“^). 
Next, using the above total errors in Great Circle coor¬ 
dinates, we converted to a and 5 through the following 
equations, 

s = — sin(iNCL) sin(Nu) sin(MU — node) 

-I-cos(incl) cos(nu), 

c = — sin(iNCL) cos(mu — node), (4) 

<7«sdss = + (5) 

and 

C^isDss = ^(s-cr^)2 + (c-cr^)2, (6) 

where mu, nu, node, and INCL refer to fields in the 
SDSS CasJobs “Frame” table. MU and NU refer to the 
Great Circle coordinates of the frame center, and iNCL 
and NODE are the inclination and right ascension of the 
ascending node of the scan Great Gircle with respect to 
the J2000 celestial equator®. 

The 2MASS PSC uses the Tycho-2 catalog to recon¬ 
struct its coordinates in the ICRS®, with accuracies be¬ 
tween 70-120 mas. 2MASS astrometric errors are re¬ 
ported as an error ellipse with the entries err_MAJ 
(o^maj), err_min (ctmin), and err_ang (erg). These 
were converted to da and as components using the fol¬ 
lowing equations, 

0’a2MAss = \/(ctmaj • sinCTe)2 -I- (ctmin • cosaeY (7) 
and 

0'52mass = Vi^MAj ■ COS ag )"^ + (ctmin ' sincre)^. (8) 

Relative astrometric calculations for 2MASS and WISE 
are computed on the unit sphere, and therefore the a 
component of the astrometric uncertainties already ac¬ 
counts for the cos S term (Vandana Desai, personal com¬ 
munication) . 

WISE is tied to the ICRS through 2MASS. However, 
to address possible systematic proper motion shifts be¬ 
tween the two catalogs due to their different epochs, the 
WISE pipeline used proper motion data from UCAC4 
to readjust 2MASS positions before they were used as 
reference stars^. Estimated errors for the source catalog 
are < 100 mas. WISE astrometric errors are denoted by 
the entries SIGRA (ctsigra) and SIGDEC (ctsigdec)- 

® https://www.sdss3.org/dr10/algorithms/astrometry.php 

® http://www.ipac.caltech.edu/2mass/r6leases/allsky/doc/ 
sec2_2.html 

‘ http://wise2.ipac.caltech.edu/docs/release/allwise/ 
expsup/sec2_5.html 


For sources with only two observation epochs, the un¬ 
certainties for each proper motion component are given 
as 


and 


/ cos <5 


V A, 


A + l^l K+<) 


-I- ct| tan^ S 


( 9 ) 




-2 




a; 


K+O 


( 10 ) 


where Aq, = a 2 — ai, A^ = 62 — Si, At = t 2 — ti, S is 
the weighted-mean declination, and ctj is the error in the 
weighted mean declination. Here, a, S and t refer to the 
position and time, with “1” being the first epoch and 
“2” the second epoch. The final term in is orders 
of magnitude smaller than 1 mas yr“^ and is therefore 
negligible. Proper motion errors for the a component are 
in proper units (i.e. Aq • cosJ). 

On average, the temporal uncertainty is much smaller 
than the time baseline between measurements {at ~ 60 
seconds for 2MASS and SDSS). The AllWISE catalog 
combines observations from the initial WISE mission 
and the post-cryogenic NEOWISE survey. We define 
the temporal uncertainty to be halfway between the dif¬ 
ference of the most recent observation and the earliest 
observation in the VFl-band (i.e. = [wImjdmax - 

w1mjdmin]/2). As is shown in FigurelH this causes some 
observations to have temporal uncertainties between 80- 
200 days. For sources with time baselines of at least four 
years, this uncertainty is 1 mas yr“^. However, a small 
fraction of our sources have time baselines of about one 
year. For this reason, we do not remove the temporal 
uncertainty term of the proper motion error. Our mo¬ 
tivation for using the AllWISE catalog rather than the 
All-Sky Release Gatalog is that the AllWISE catalog has 
better astrometric accuracy. This is due to the inclusion 
of proper motions to correct 2MASS astrometric refer¬ 
ence stars for the greater than nine year baseline between 
the WISE and 2MASS surveys. 

For sources with three observational epochs (i.e. 
WISE, SDSS, and 2MASS), we computed a weighted lin¬ 
ear ht to the positions versus time. Rather than use an 
linear least squares approach, which requires the uncer¬ 
tainty in the independent variable (in this case time), to 
be negligibl e, we chose to invoke an Orthogonal Distance 
Regression (|Boggs & Rogers||1990| method to calculate 
proper motions in each component {a, S) separately. This 
allows us to take into account the sometimes significant 
temporal uncertainties on WISE observations. 


2.2. Precision: Measuring the Motions of SDSS 
Quasars 

To investigate the intrinsic error in the proper motions 
for each survey, we required stable objects on the sky 
with essentially zero tangential motion. Quasars make 
ideal calibrators due to their high luminosities and ex- 
tragalac tic distances. We used the DRIO SDSS quasar 
catalog ( Paris et'ar]|2014 ), which contains 166,583 spec- 
troscopically conhrined quasars. Using this catalog, we 
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Figure 1. WISE temporal uncertainties for our sources. The 
large uncertainties are due to the AllWISE catalog being a com¬ 
posite of WISE and NEOWISE observations. This temporal un¬ 
certainty will account for < 1 mas yr~^ in our astrometric solution 
for stars with baselines > 4 years, however, we account for it since 
a small subset of our observations have baselines of ~1 year. 


Table 1 

Limits for WISE, SDSS, and 2MASS 


Figure 2. Unique SDSS QSO matches per search radius in 
2MASS and WISE. Both reach a maximum at a search radius of 
0.5'L The increasing number counts at > 2" for 2MASS indicate 
resolved neighboring objects being pulled into our search radius for 
faint QSOs that 2MASS is not able to detect. 

9. The closest neighboring primary SDSS object was 
greater than 6.1" (the W1 FWHM) from our 
source. 


Survey 

Band 

Limiting Magnitude 
(mags) 

Limiting Flux 
(ergs s“^ cm“^ A 

SDSS 

r 

22.2 (AB) 

~4 X 10-1® 

SDSS 

i 

21.3 (AB) 

~6 X 10-1® 

SDSS 

z 

20.5 (AB) 

~8 X 10-1® 

2MASS 

j 

15.8 (Vega) 

~2 X 10-11 

2MASS 

H 

15.1 (Vega) 

~4 X 10-11 

2MASS 

Ks 

14.1 (Vega) 

~8 X 10-11 

WISE 

W1 

17.1 (Vega) 

~io-i® 


cross-matched to both the 2MASS PSC and WISE All- 
WISE source catalog. Matching was done using search 
radii in steps of 0.5" out to 4", the results of which 
are shown in Figure The number of unique matches 
reaches a maximum at 0.5" for both 2MASS and WISE 
matches. Since we sought to estimate the precision of our 
catalog, we chose to only use matches within a search ra¬ 
dius of 0.5". This left us with 69,949 matches between 
SDSS and WISE and 2351 matches between SDSS and 
2MASS, 2283 of which were common in both WISE and 
2MASS. The low number of matches between SDSS and 
2MASS is due to the difference in wavelengths and rela¬ 
tive depths between the two surveys. The mamitude and 
flux limits of each survey are shown in Table Because 
we were assessing the intrinsic precision of the catalog, 
we chose to apply further cuts to retain only the most 
pristine detections of QSOs for our analysis. The follow¬ 
ing criteria were required: 


These cuts left us with 447 matches for sources with 
1F/5F+SDSS-I-2MASS, 4091 matches for sources with 
only WISE+SDSS, and 413 matches for sources with 
only SDSS-h2MASS. 

We applied the above algorithm to compute the angu¬ 
lar distance measured from the surveys. Distributions 
of our computed angular distances are shown in Fig¬ 
ure]^ The estimated errors in our sample of QSOs are 
in agreement with, or slightly better than, reported po¬ 
sitional uncertainties among each of the three surveys. 
The largest uncertainties are in our SDSS-I-2MASS base¬ 
line, however, these objects make up a small fraction of 
our entire proper motion catalog (see Section §. Fig¬ 
ure ^ represents the intrinsic positional errors in each of 
our fits (~90 mas for FF/S'A-I-SDSS-I-2MASS, ~80 mas 
for WISE+SDSS, and -125 mas for SDSS-H2MASS). We 
add this error, weighted by the time baseline, to the 
proper motion error for each component. 

A number of studies have attempted to correct system- 
atic errors in proper motion catalogs using QSOs (e.g., 


Roser et al. 120101 Wu et al.|2011[ Lopez-CorredoirapOld 

Grabowski et al.|2015D. We did not observe any large sys- 

tematics within the WiSE+SDSS sources, and the small 
sample size among all other surveys did not allow us to 
investigate or correct for any potential systematics off¬ 
sets. 


3. DATA 


1. SDSS CLEAN = 1, 

2. 2MASS CC_FLG = 000, 

3. WISE CC_FLAG = 0000, 

4. WISE wlSNR >= 30, 

5. 2MASS J_SNR, H_SNR, or k_snr >= 10, 

6. 2MASS GAL.CONTAM = 0, 

7. WISE EXT_FLAG = 0, 

8. tjQ < 175 mas & cr^ < 175 mas, and 


3.1. Building the Catalog 

SDSS DRIO boasts over 900 million unique optical 
sources (after accounting for multiple epochs of photom¬ 
etry)®. This study focused on the low-mass stars (point 
sources) within DRIO. We used the following photomet¬ 
ric selection criteria, many of which were adapted from 


https://www.sdss3.org/drlO/scope.php 
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WISE + SDSS + 2MASS WISE + SDSS SDSS + 2MASS 

Aa cos 6 (mas) Aa cos 6 (mas) Aa cos 6 (mas) 



Figure 3. Normalized distributions for the positional difference components among each of the surveys (gray). The best-fit normal 
distribution is also plotted (red dotted line). The precision tends to be better for observations with a WISE epo ch. These uncertainties 
are weighted by the time baseline between observations and added in quadrature to the fitting errors (Section|2.1[). 


Bochanski et al. 
search: 


(2010 


hereafter BIO), for our initial 


1. The photometric objects were flagged as primary. 
This was done by querying the Star sub-catalog 
of the PhotoPrimary catalog. This also ensured 
a morphological classification of a point source. 

2. The photometric objects fell within the following 
color limits: 

16 < r < 22, 
i < 22, 
z < 21.2, 
r — i ^ 0.3, and 
i — z ^ 0.2. 

These color criteria allow for sources slightly bluer 
than typical dM colors, but were chosen to be in¬ 
clusive for this stage of the selection process. The 
i- and z-band magnitude limits extend past the 
SDSS 95% completeness {i < 21.3 and z < 20.5) 
limits, but we apply additional cuts below. The r- 
band limit removes saturated sources and is slightly 
brighter than the 95% completeness limit (r < 
22.2). Although we obtain sources past the 95% 
completeness limit, we apply more stringent crite¬ 
ria later. 

3. We removed sources that had a flag in the r-, i-, or 
z-bands indicating: 


(a) saturated photometry (saturated); 

(b) a signihcant amount (>20%) of the flux was 
interpolated from the point spread function 
(PSF; psf_flux_interpreted) 

(c) centroiding failure caused center to be deter¬ 
mined by peak pixel (peakcenter); 

(d) too few good pixels for an interpolated 
source, causing errors to be underestimated 
(bad_counts_error); 

(e) center pixel was too close to interpolated pixel 
(iNTERP. center), and included pixels that 
were not checked for peaks, potentially sat¬ 
urated (notchecked); 

(f) after deblending the object did not have a 
peak (deblend_nopeak); and 

(g) contained a pixel interpreted to be part of a 
cosmic ray (cosmic_RAY). 

4. Sources that indicated centroid was not determined 
from the r-band, but were transformed from some 
other band (canonical.CENTEr), were removed. 
This was done to ensure objects with reliable r- 
band astrometry, which was used as the calibrator 
for our SDSS baseline. 

These cuts ensured that the riz photometry for each 
source were reliable. Applying these criteria returned 
69,792,454 objects. 

Using a 6 " search radius (a similar mat ching radius 
to other searches for ultra-cool dwarfs, e.g., [Zhang et al. 
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20091, we matched the SDSS photometric stars to WISE 
and 2MASS sources requiring: 


1. 2MASS J_PSFCHI, H_PSFCHI, Or K.PSFCHI < 2. 
This ensures a point-like source morphology, and 
increases our confidence for a good astrometric 
measurement. 


2. WISE w1rchi2 < 3. This is the same require¬ 
ment used to determine a single-point source from 
a blended object (after deblending) by the WISE 
photometric pipeline. 

3. WISE SNRvyi > 0. This removes sources that 
were not detected in Wl, but detected in a longer 
wavelength band. The spectral energy distribu¬ 
tions (SEDs) of low-mass stars peak in the near- 
infrared, therefore, as Wl was the deepest of the 
WISE bands, low-mass stars should have a detec¬ 
tion in W1 if they have a detection at a longer 
wavelength. 

After this cut we were left with 20,164,221 matches 
with entries in all three catalogs, 22,741,703 with 
only 14757?-I-SDSS detections, and 2,947,606 with only 
SDSS-I-2MASS detections. 


3.2. Tracing the Stellar Locus 

We were interested in selecting point sources that have 
the expected colors of main-sequence stars for our initial 
sample. To ensure this, we computed the expected main- 
sequence colors for low-mass and very-low-mass stars in 
the WISE, SDSS, and 2MASS photometric s ystems. We 
requir e d the following criteria, adapted from 


|Davenport| 
IWISE Ex- 


et al. (2014) for low-mass stars, and the All 
planatory Supplement®, to select high quality matches: 


1. 13.8 < r < 21.5 & < 0.05, this cut ensured 

precision photometry within the SDSS 95% com¬ 
pleteness limit and saturation limit, and ensured 
good morphological classification, which has been 
shown to have an error rate of 5% at r = 21 (Lup- 
ton et al.||2001 |. 


2. J > 12 & aj < 0.05, this cut ensured precision 
photometry, and should have removed gian t stars 
by using a lower limit on the magnitude (Covey 
et al.|[2^ . 


3. W1 < 17.1 & awi < 0.05, this cut also ensured 
precision photometry and selected only sources 
within the AllWISE 95% completeness limit. 


4. |6| > 20°, this cut reduced extinction effects by 
removing candidates near the Galactic plane. 

After applying the aforementioned cuts, we were left with 
9,298,344 sources with only SDSS photometry, 9,890,521 
sources with 1475E-I-SDSS photometry, and 2,390,962 
sources with SDSS-I-2MASS photometry. 

Previous studies have suggested that the majority of 
point-like sources that meet our color selection criteria 


^ http://wise2.ipac.caltech.edu/docs/release/allwise/ 
expsup/sec2_4a.html 


(Section |3.1|) a 
ski et al^ 2010 


(Section m. are stars versus distant Galaxies (|Bochan?| 


10) , however, we wish to be more selective, 


choosing only high probability stellar candidates for our 
catalog. To do this, we took our above candidates, se¬ 
lected for their reliable photometry, and examined the 
stellar locus in bands from all three surveys. Many pre¬ 
vious studies have investigated the stellar color locus in 
num erous colors and a num ber of photometric systems 
(see Davenport et al. 2014 and references therein). In- 
stead of using previous results, we chose to measure the 
color locus from our stellar sample chosen above. We 
also chose to only measure the locus for the reddest SDSS 
bands (riz) and the deepest 2MASS and WISE bands (J 
and Wl, respectively). Our chosen colors for computing 
the stellar locus were r — i versus i — z, r — z versus z — J, 
and r — z versus z — Wl. 

In small steps of (5(color), we computed the median 
absolute deviation from the median and removed objects 
greater than 5 times the median of the deviations until 
our distribution converged (i.e. there were no further 
sources to remove). We then computed the mean and 
standard deviation of the remaining color distribution. 
We used a bin size of 0.01 mags for the high density areas 
(the middle of the distribution), taking steps of 0.01 mags 
and recomputing. For the red end of the distributions 
(the low density areas), we increased the bin size to 0.2 
mags. Averages and Icr colors, color steps, bin sizes, 
and the number of stars in each bin can be found in 
Appendix ^ The visual representation of our computed 
means ancTScr deviations are shown in Figure]^ where it 
can be seen that we trace the source density, removing 
large color outliers. 

Due to our initial selection criterion for sources along 
low-extinction sight-lines (|&| > 20°), it is possible that 
our sample could be biased towards stars with low extinc¬ 
tion values. Rather than apply an extinction correction 
to our stars, we investigated how extinction may bias our 
selection criteria. To test how accurate our color criteria 
are for selecting low-mass stars and the effects of inter¬ 
stellar reddening, we applied our cuts to th e SDSS DR7 


spect roscopic sample of 70,841 M dwarfs (West et al. 
201 1| hereafter Wl l), all of wh ich have estimated ex- 


tmction values from Jones et al. (2011). We investigated 
the fraction of returned Wll stars, after passing them 
through our color selection crite ria. We applied the cuts 
above and the cuts from Section ^ to the Wll catalog, 
excluding the requirement that stars have |&| > 20° since 
we were interested to see how extinction affected our se¬ 
lection method. The total return fraction for stars with 
Av < 0.5 and Ay > 0.5 are shown in Table for each 
combination of color selection criteria. For stars with 
Ay < 0.5, our color selection criteria returned > 95% 
of the Wll inputs stars for all color selection criteria. 
Even for stars with Ay > 0.5 our color selection criteria 
returned more than 92% of the Wll input stars. 

The Wl l catalog also contai ns both disk dwarfs and 
subdwarfs (jSavcheva et al.|2014|), populations that can be 
separated b y their metalljcity u sing the ^-spectroscopic 
index (e.g., Dhital et al. 2012). To test if our catalog 
preferentially selects disk dwarfs (C > 0.825) or subd¬ 
warfs (C < 0.825), we again investigated the fraction of 
returned Wll stars after passing them through our color 
selection criteria, and results are shown in Table 
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Table 2 

Color Selection Criteria Returns for Wll 


Color Criteria 

Total 

Disk Dwarfs 

Subdwarfs 


Av 

< 0.5 


la 

98.8% 

98A% 

994% 


99.2% 

99.2% 

98.7% 

3'^ 

99.6% 

99.6% 

95.8% 

l-b2 

98.3% 

98.3% 

98.3% 

2-1-3 

99.0% 

99.0% 

95.0% 

l-b3 

98.5% 

98.6% 

95.8% 

1+2+3 

98.1% 

98.1% 

95.0% 


>lv 

> 0.5 


la 

98.4% 

984% 

984% 

2*> 

97.2% 

97.2% 

96.3% 

3'^ 

97.0% 

97.2% 

94.4% 

1+2 

96.2% 

96.2% 

96.3% 

2+3 

95.9% 

96.2% 

92.6% 

1+3 

96.1% 

96.2% 

94.4% 

1+2+3 

94.9% 

95.2% 

92.6% 


— ivs. i — z criteria. 

^ r — z vs. z — J criteria. 
r — z vs. z — W1 criteria. 


Our selection criteria are well-suited for retrieving both 
disk dwarf and subdwarf populations, with selection of 
stars in high extinction environments being slightly less 
reliable, but all above 92%. Our ability to select the vast 
majority of spectroscopic low-mass stars based on colors 
alone demonstrate that these color criteria are suitable 
for selecting photometric low-mass stars, even in regions 
of moderate extinction. 

3.3. Initial Stellar Sample 

Applying the color cuts from Figurej^to our initial stel¬ 
lar sample, we were left with 24,571,934 stars. The color 
and magnitude range over which each of our matched 
samples is found is shown in Figure lF75i?-|-SDSS 
sources are typically fainter and redder, biasing these 
matches to later-type stars. SDSS-I-2MASS matches are 
bluer, biasing these matches to the earliest type stars in 
our sample. 

Using the above sources as our input cata log, we com¬ 
puted proper motions following Section [O] The surveys 
that went into each proper motion measurement (e.g., 
WT5A-I-SDSS-I-2MASS or WTS'F-I-SDSS) are indicated 
in our catalog by a three bit flag (dbit), where each bit 
represents if a survey was used in the proper motion mea¬ 
surement. Descriptions of dbit and the number of initial 
sources with each fit are shown in Table |3l 

3.4. Reducing Contamination 

There are two main issues that can contribute to spu¬ 
rious motion estimates for our stars: 1) nearby sources 
that can offset the center of the measured PSF if de¬ 
blending methods fail; and 2) short time baselines that 
can contain motion caused by parallactic effects in addi¬ 
tion to tangential motion. We will closely examine both 
of these issues in the following subsections. 

3.4.1. Nearby Neighbors 

Close neighboring sources are expected to be a con¬ 
taminant, especially for WISE stars due to the large 
beam size (FWHMi4/i=6.1"). WISE active deblending 
allows, at most, two components to the PSF fit. The 





Figure 4. Density plots in WISE, SDSS, and 2MASS colors, each 
bin is (0.01 mag)^. The first filled contour is drawn at 10 sources 
per bin, and each subsequent filled contour increases by a factor 
of 10. Mean colors and 3 (t errors used to select our initial stellar 
sample are also shown (dotted line and solid lines, respectively). 
Our method well approximates the stellar locus. The points used to 
draw the limits (dashed and dotted lines) are shown in Appendix[A| 
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Figure 5. Hess diagram (r vs. r — z) with density contours, 
each bin is (0.1 mag)^. W75'S+SDSS+2MASS are drawn in blue, 
WISE-\-ST)SS are drawn in red, and SDSS+2MASS are drawn in 
green. Contours are drawn at 1000, 5000, and 10000 stars per bin 
(solid, dash-dotted, and dashed lines, respectively). Each subsam¬ 
ple is biased towards a specific brightness and color range. 


Table 3 

Proper Motion Detection Flags 


DBIT 

Description 

Number 

on 

SDSS and WISE surveys were used 
to calculate the proper motions. 

initial'. 11,911,109 
final: 1,801,369 

no 

2MASS and SDSS surveys were used 
to calculate the proper motions. 

initial: 1,052,228 
final: 69,199 

111 

All three surveys were available and 
used in the linear fit to calculate 
proper motions. 

initial: 11,608,597 
final: 1,733,104 

000 

SDSSH-USNO-B proper motions 
are available. 

final: 6,620,838 
unique: 5,216,855 


Note. — We show both the initial number of sources that met our first 
selection criteria (see Section and the final number of stars that went 
into our catalog. 


robustness of this deblending is likely dependent on the 
flux difference between the two blended sources. To in¬ 
vestigate how neighboring objects affect proper motion 
measurements, we used SDSS CasJobs to select all neigh¬ 
boring primary objects within 15" of our sources’ SDSS 
positions. Next, we selected only sources that had one 
neighboring object within the 15" search radius. We ex¬ 
pect sources whose photocenters have been moved signif¬ 
icantly in WISE will have larger fitting^errors for sources 
detected in all three surveys. Figure q shows fitting er¬ 
rors as a function of distance to a neighboring object and 
r-band magnitude difference between our source and the 
neighboring object. Fitting errors: (1) are significantly 
larger for brighter neighboring objects, and objects at a 
distance < 8"; and (2) decrease for extremely close (and 
bright) neighboring objects (< 2") due to a reduction in 
the offset of the measured photocenter from its true po¬ 
sition, making the measured positional deviation small. 

Figure illustrates that sources with neighboring, 
bright objects have problematic motion estimates that 
can be seen in the fitting errors. The following ques¬ 
tions must be answered to determine which proper mo- 


Figure 6. Total proper motion fitting errors for the 
W75'i?-|-SDSS+2MASS proper motions as a function of distance 
to nearest SDSS primary object and magnitude difference between 
our source and its closest neighbor. Only stars with one neigh¬ 
bor within the search radius are shown, each bin is 0.2 mags X 
0.2". Stars with fitting errors > 40 mas yr“^ are suspect due 
to blending. Dashed and dotted lines correspond to criteria, be¬ 
low which we recomputed proper motions or removed sources (see 
text for further details). For VF/SF-I-SDSS-I-2MASS stars, we re¬ 
computed proper motions for stars within the dashed lines. For 
WTSF-I-SDSS sources, we removed stars that had a neighboring 
object within < 8" (dotted line). 


tion measurements are reliable: 1) at what distance does 
a neighboring object affect our measurement; 2) at what 
magnitude difference does a neighboring object affect our 
measurement; and 3) what is the fitting error threshold 
above which a measurement is considered affected (by 
either neighbors or parallactic effects)? We will explore 
each of these three questions below. 

We expect the error distributions for sources with a 
neighboring object at a distance < 8" and > 8" to be dif¬ 
ferent due to the inclusion of more large errors for the dis¬ 
tribution with a neighboring object < 8". We require a 
statistical argument in choosing at what distance a neigh¬ 
boring object affects our measurement. To test for simi¬ 
larity (or difference) in the error distributions, we chose 
to use the Anderson-Darling test since it is more sensitive 
to the tails of the distribution, where we expect the larger 
fitting errors to reside. To determine at which distance 
a neighboring object affects our proper motion measure¬ 
ment, we selected all sources with r-gource—''’neighbor > —2. 
The requirement for sources with rgource — ^neighbor > — 2 
is relatively arbitrary, we are selecting sources where we 
know fitting errors begin to increase, later we make a 
more rigorous estimate for rgource — ''’neighbor- Next, we 
binned fitting errors in steps of 0.5" between 4" to 12", 
performing an Anderson-Darling test between adjacent 
0.5" bins, searching for neighboring distributions with 
the largest dissimilarity (as traced through a minimum 
in the p-value). We found a minimum p-value (8 x 10“®) 
between the bins from 7.5" to 8" and 8" to 8.5", making 
8" our cutoff value. 

Next, to determine at which magnitude difference 
(''’source—^neighbor) a neighboring object affects our proper 
motion measurement, we selected all sources with a 
neighboring object < 8". Again, we binned fitting errors 














The MoVeRS Catalog 


9 


in steps of 0.5 mag between —5 to 3 mags, and performed 
Anderson-Darling tests between adjacent distributions. 
We found a minimum p-value (4 x 10“^) between the 
bins from Csource — ^'neighbor = ”3 to —2.5 and —2.5 to 
—2, making —2.5 our cutoff value. Additionally, we vi¬ 
sually inspected a number of stars with proper motion 
errors between = 30-50 mas yr“^ and found that 
sources with > 40 mas yr“^ tended to be spurious 
measurements. We chose to recompute proper motions 
for these stars using just two of the three surveys. De¬ 
tails regarding which surveys were chosen ( WISE+SDSS 
or SDSS-I-2MASS) are described below. 

Our method for recomputing proper motions is as fol¬ 
lows: 

1. We completely removed stars with a neighboring 
source found within 8" and with rsource~*’neighbor > 
-2.5. 


2. For stars that do not meet the above criteria, but 
which have > 40 mas yr“^ in either the a or S 
component, we recomputed proper motions by first 
requiring a minimum signal-to-noise (S/N) in each 
survey. For WISE, we required S/N^i > 3, and for 
2MASS, we required S/Nj^h, or a, > 5. If only one 
survey met our minimum S/N threshold, we use 
that survey. If both surveys met our S/N threshold, 
we used the survey with the highest S/N. To use 
the WISE baseline, we further required a minimum 
distance to the nearest neighbor be > 8". 

(a) If all the previous criteria of (2) were met, 
and the highest S/N between both surveys 
was equal, we used the two surveys with the 
longest time baseline (either W/S'A-I-SDSS or 
SDSS-k2MASS). 

3. Lastly, we required a time baseline > 1 year be¬ 
tween observations. Time baselines shorter than 
this are susceptible to parallax effects for nearby 
stars (oo20 mas yr“^ for a star at 100 pc with a 
baseline of 6 months). 

As a final note, the buildup of large errors at a neigh¬ 
boring distance of may also be due to our search ra¬ 
dius (6") picking up neighboring objects as the primary if 
our source becomes too faint at longer wavelengths. Re¬ 
moving these sources ensures both possibilities for con¬ 
tamination are removed from our catalog. Since we are 
only interested in stars exhibiting bona-fide tangential 
motions, we kept only stars with total proper motions 
greater than twice their uncertainty (p.tot > 


4. RELIABILITY OF PROPER MOTIONS 

To assess the reliability of our catalog, we com¬ 
pared our proper motio ns to the LSPM catalog , the 
SDSS-I-USNO-B catalog (|Munn et al.||2004| |2008 here- 
afer M04), and the recent de ep survey complet ed within 
a 1098 deg^ SDSS footprint ( Munn et al.|2014 hereafter 
M14). The reliability of LSRAI should be close to 100% 
as all these sources have been verified by eye. LSPM 
stars are selected for larger proper motions (> 150 mas 
yr“^), and the catalog is biased towards brighter stars 
due to the use of Schmidt plates for the earliest base¬ 
lines. M04 has time baselines of ~50 years, giving it a 


high precision, but M04 is not as deep as our catalog, 
also due to the use of Schmidt plates. M04 has also been 
matched to a high number of SDSS sources since SDSS 
was used as the most recent baseline in computing proper 
motions. Lastly, M14 allows us to test the fidelity of our 
faintest sources, which LSPM and M04 do not probe. 
Together, these catalogs allow us to assess the reliability 
of our sources for the stars with both small and large 
proper motions, and across all magnitudes. 


4.1. Comparison to LSPM 


The LSPM catalog contains 61,977 stars in the north¬ 
ern hemisphere with proper motions > 150 mas yr“^. 
The precision of LSPM is ~8 mas yr“^. LSPM was 
not specifically designed to target low-mass stars, how¬ 
ever, due to its selection of high proper motion stars, it 
primarily consists of nearby dwarf stars. We matched 
our catalog to LSPM stars using their 2MASS desig¬ 
nations, which produced 12,930 matches. We investi¬ 
gated the agreement between our catalog and LSPM in 
all three subsamples (e.g., 1T/S'A-|-SDSS-I-2MASS) for 
each proper motion component, as shown in Figure!^ A 
small number of large outliers were identified between our 
1T/S'i?-|-SDSS-|-2MASS and LSPM matches (Figure Jj^ 
red circles). We looked through archived images (D^, 
SDSS, 2MASS, and WISE), since these stars should all 
have apparent proper motions over the > 20 year base¬ 
line in archived images. The majority of these stars 
showed proper motions more consistent with our mea¬ 
surements in the archived images. The cause of the spuri¬ 
ous measurement in LSPM is unclear, presumably these 
are bad residual images in the SUPERBLINK pipeline 
that eluded inspection. Many of these stars also had 
proper motions consistent with our measurement in an¬ 
other catalog (e.g., USNO-B, M04, NLTT, PPMXL, or 
URATl). 

We compared the 2cr agreement in both proper mo¬ 
tion components to test the absolute agreement between 
our catalog and LSPM. For our 1T/S'i?-|-SDSS-|-2MASS 
matches, we found 98% agreement between our catalog 
and LSPM. For our WTS'A-I-SDSS and SDSS-I-2MASS 


matches this agreement was 97% and 96%, respec¬ 
tively. All agreements increase to 99% at 3tT. All our 
WISE-\-SD SS ma tches to LSPM were recomputed values 
(see Section 3.4.1), since these stars are all bright enough 
to have an entry in 2MASS. This comparison shows our 
catalog is reliable for the fastest moving stars. 


4.2. Comparison to MOf 

Due to the sensitivity of M04, spurious proper motion 
estimates are expected for fainter stars. To compare our 
catalog to M04 we needed to determine w hich sources 
were r eliable. We matched our catalog to the |Munn et al.| 
(20041 catalog (contained within SDSS CasJobs as the 
ProperMotions table) by SDSS objid. We chose to 
use only sources wit h “good” proper m otions, adopting 
the constraints from Kilic et al. (2006): 


1. dist 22 > 7, the nearest neighboring objects with 
g < 22 is more than 7" away. 


2. MATCH = 1, there is a one-to-one match with the 
USNO-B object and the SDSS object. 
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1200 




Figure 7. Residual proper motion (/ithis study ~ Mlspm) for both components (/^.q; and fis) as a function of our total proper motion. 
Typical errors are shown in the top right corners. Top: 1T/5F+SDSS+2MASS stars compared against LSPM. Overall, our stars show 
good agreement with LSPM, however, a small number of outliers were identified and investigated (red circles, see text for details). Middle: 
WISE-\-SDSS stars compared to LSPM. Again, there is good agreement between our stars and LSPM, with the outliers typically due to 
bad astrometry in one of our surveys. Bottom: SDSS+2MASS stars compared to LSPM. Similar to our WISE+SDSS stars there is good 
agreement with LSPM, with the outliers once again due to bad astrometry in one of our surveys. 
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3. nFit = 6, the object was detected in all five USNO- 
B plates, as well as detected in SDSS 

4. SIGRA < 525, and SIGDEC < 525, the RMS resid¬ 
uals for the proper motion fits were less than 525 
mas in both components. 


These criteria yielded 919,867 matches between our 
SDSS sources and M04 sources. 

Since M04 is biased towards brighter objects, 
we chose to compare it to sources that have 
proper motions measured using all three surveys 
(lT/5'£’-|-SDSS-l-2MASS). There were 842,776 matched 
stars with 1T/6'E-|-SDSS-|-2MASS measured proper mo¬ 
tions, making up the bulk of the matches between our 
catalog and M04. We investigated the residuals between 
our computed proper motions and those from M04 as 
a function of color and magnitude (Figure |^. There 
is a strong increase in the magnitude of the residuals 
for fainter stars and bluer stars. The correlation be¬ 
tween apparent magnitude and r — z color is due to 
bluer stars peaking in the optical (and NIR), while get¬ 
ting fainter at WISE bands. The combined filter set 
of 1F/S'F-|-SDSS-|-2MASS is more sensitive to stars that 
peak towards the red end of the NIR, out to ^ 4 /rm. We 
expect M04 to be more sensitive to bluer stars, however, 
we can examine the proper motion distributions for these 
faint sources to determine whether ours, or the M04 val¬ 
ues are more consistent with expectations. 

In general, we expect fainters stars to be more dis¬ 
tant, and therefore exhibit smaller proper motions. Our 
proper motions for fainter sources tend to be higher than 
those from M04, as can be see in Figure (positive 
residuals). We performed a linear fit to the residuals 
in color-magnitude space, fitting to where the residu¬ 
als became higher than 20 mas yr“^(dashed line. Fig¬ 
ure]^, and where the residuals became higher than 10 
mas yr“^(dash-dotted line. Figure]^. Our results are 
shown in Figure with the proper motions separated 
by, 

r < 16.65-1-1.435(r—z) for A/rtot < 20 mas yr“^, (11) 
and 

r < 15.34-1-1.896(r—z) for A^tot < 10 mas yr“^. (12) 


As stated above, due to the difference in wavelengths 
between each of the surveys used in calibrating our cat¬ 
alog, we expect fainter, redder sources to be more eas¬ 
ily detectable than fainter, bluer sources. This is ex¬ 
hibited in Figure where our proper motion errors also 
begin to grow large at approximately the same limit that 
our residuals began to increase. Therefore, both our 
proper motions and our proper motion errors increase at 
the same limit, potentially keeping our measurements to 
within errors but with less precision. To quantify the reli¬ 
ability of our catalog, we compared the 2a and 3 ct agree¬ 
ment for both proper motion components for ma tched 
stars, both fainter and brighter than equation (ITR, and 
as a function of color, the results of which are shown in 
Table H 

In general, there is good agreement between our mea¬ 
surements and the M04 catalog. This agreement is less 
certain for stars fainter than equation (11), but increases 



Figure 8. Top: Residual total proper motion (A/rtot = 
Ttot (wrSB+SDSS+2MASS) “ Mtot {M04) 3'® ^ function of source r- 
band magnitude and r — z color, each bin is (0.1 mags)^. The 
performance of our catalog is tightly correlated with source color 
and magnitude. The dashed line represents the limit above which 
residuals become > 20 mas yr“^, and the dash-dotted line rep¬ 
resents the limit above which residuals become > 10 mas yr~^. 
More stringent criteria used to investigate the reliability of the 
M04 catalog (Section]^ are also shown (dotted lines). Bottom: 
Total proper motion error as a function of source r-band 

magnitude and r — z color, each bin is (0.1 mags)^. The dashed 
and dash-dotted lines are the same as the top. Our measurement 
uncertainty becomes large at the same limit that residuals begin 
to grow larger than 20 mas yr“^. The performance of our catalog 
is tightly correlated with source color and magnitude, due to our 
inability to probe small proper motions expected for more distant 
(fainter) sources. Although we lose precision at fainter magnitudes, 
our proper motion measurements are still consistent to within the 
errors. 


for redder colors fainter than this limit. Since the per¬ 
formance of the USNO-B survey is also expected to de¬ 


teriorate at fainter magn itudes (V > 19; Lepine & Shara 


2005 Dong et al. 2011), we can expect sources in dis- 


agreeinent to be a mixture of spurious measurements in 
our catalog and in USNO-B. This is left as a caution¬ 
ary note, as many of the stars fainter than equation (111 
will have real proper motions. We do not remove these 
stars, as many of them will be true measurements, we 
instead augment them with a second measurement from 
M04 (see Section]^. 
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Table 4 

Agreement with other Proper Motion Catalogs 


Comparison Catalog 

Subsample 

r — z 

2a Agreement 

3a Agreement 

M04 

WISE+S'DSS+2MASS 

All 

86% 

95 % 

M04 

SDSS-I-2MASS 

All 

28% 

72% 

M04 

WISE+SDSS 

All 

15% 

50% 

M14 

WISE+SBSS+2MASS 

All 

82% 

94% 

M14 

SDSS-I-2MASS 

All 

21% 

74% 

M14 

WISE+SDSS 

All 

25% 

63% 

M14 

WISE+SDSS 

> 3 

82% 

93% 


r < 

16.65 -1- 1.435(r - z) 



M04 

WISE+S'DSS+2MASS 

All 

94% 

99% 

M04 

SDSS-I-2MASS 

All 

87% 

97% 

M04 

WISE+SDSS 

All 

46% 

63% 

M14 

WISE+SDSS+2MASS 

All 

95% 

99% 

M14 

SDSS-I-2MASS 

All 

69% 

95% 

M14 

WISE+SDSS 

All 

57% 

77% 

r > 16.65 + 1. 435 ( 7 * — z) 

M04 

WISE+SDSS+2MASS 

All 

46% 

75% 

M04 

WISE+SDSS+2MASS 

> 1 

50% 

77% 

M04 

WISE+SDSS+2MASS 

> 1.5 

52% 

78% 

M04 

WISE+SDSS+2MASS 

> 2 

55% 

80% 

M04 

SDSS-I-2MASS 

All 

13% 

66% 

M04 

WISE+SDSS 

All 

13% 

49% 

M04 

WISE+SDSS 

> 2 

28% 

59% 

M14 

WISE+SDSS+2MASS 

All 

60% 

84% 

M14 

WISE+SDSS+2MASS 

> 3 

93% 

98% 

M14 

WISE+SDSS 

> 3 

80% 

92% 


A similar comparison can be made between our 
H/T5i?+SDSS and SDSS+2MASS proper motions and 
M04. However, these samples are limited in that 
they are typically bluer (SDSS+2MASS) or fainter 
(WTS'A’+SDSS) than our 1F/S'i?+SDSS+2MASS stars. 
These comparisons are also shown in Table In gen¬ 
eral, agreement is again better for stars bri^ter than 
equation ([II| ), and for redder stars. 

We further investigated how our agreement scaled with 
proper motion errors among all three subsamples, as is 
shown in Figurej^ Stars with larger proper motion errors 
(e.g., fainter stars and stars with smaller time baselines) 
tend to have to be closer to our Scr limit, or be in disagree¬ 
ment. The majority of stars across all three subsamples 
tend to be within Scr agreement, a proper motion error 
cut may be useful in selecting higher confidence proper 
motions. It is difficult to say if the disagreement between 
the faintest sources in our catalog (and M04) is a limi¬ 
tation of M04 or our catalog. To answer this question, 
we require a proper motion catalog as deep (or deeper) 
than our catalog. 

4.3. Comparison to MI 4 

The M14 proper motion survey covers only a small 
footprint within SDSS, but provides proper motions 
down to r « 22. This catalog has time baselines > 5 
years with a precision of ~10 mas yr“^. We applied all 
the M14 suggested quality cuts (summarized in Table 1 of 
M14). Although these cuts were designed to remove spu¬ 
rious proper motion measurements, M14 found that they 
were only able to remove ^50% of the “bad” detections, 
while retaining >97.6% of the “good” detections (from 
comparison stars between M14 and LSPM). Therefore, 
it is possible that some of the proper motions retrieved 


are spurious measurements. The hdelity of the proper 
motions within this catalog have not been independently 
verified, therefore, the following analysis should be taken 
as both an investigation into the reliability of our catalog 
as well as M14. 

We matched our catalog to M14 by SDSS OBJID, giving 
us 510,799 matches. We performed an analysis similar to 
our comparison to M04 above. We began with a compar¬ 
ison of our IF/S'i?-|-SDSS-|-2MASS stars, which made up 
246,678 of our matches. Figure [T0| shows the total proper 
motion residuals (/^tot (w/SE+Sdss-i-2 MASS) ~ Mtot (M14)) 
as a function of r-band magnitude and r — z color. Equa¬ 
tions ( [TI| ) and ( [1^ from our M04 analysis represent sim¬ 
ilar reliability for M14, past the magnitude limit that 
M04 was able to probe. Reliability tends to increase for 
stars with r — z > 3. 

We again explored the 2cr and 3cr agreement be¬ 
tween both proper motion components for M14 and 
our catalog (see Table |^. Agreements between our 
lF7S'i?-|-SDSS-|-2MASS measurements and M14 are simi¬ 
lar to those found with M04, except for the HTS'F-I-SDSS 
stars which are in better agreement than those found 
against M04 with a larger number of stars (261,141). 
Agreement is also better for redder stars (stars with 
r — z > 3), the majority of w hich were fainter than 
equation (111 from Section 4.2 This shows that our 
WISE+SDUti sources are typically too faint to be have 
reliable proper motions in M04. In general, our proper 
motions across all three subsamples appear to be reli¬ 
able (to within our uncertainties) up to our magnitude 
limit of r = 22 , w ith redder sources having higher reli¬ 
ability. Figure [D shows how proper motion agreement 
scales with proper motion error (similar to Figure |^. 
We find slightly better agreement for stars with larger 
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WISE+SBSS+2MASS 



SDSS+2MASS 



WISE+SBSS 



Figure 9. Violin plots for residual total proper motion (Aptot = Ptot {1T/SB+SDSS+2MASS) ~ Mtot (M04)) ^ function of total proper 

motion error. These plots show the relative distribution of stars for each error bin. The dashed and dotted line represent 2cr and 3 (j 
agreement, respectively. Top: The majority of IT/5'£^+SDSS+2MASS stars are all in 3cr agreement with M04, agreement is better for 
small errors. Middle: The SDSS+2MASS stars tend to be in better agreement for smaller errors. From our error distribution (Section!^, 
the majority of our sources should have reliable (within Sct) proper motions. Bottom: The WISE-\-SDSS stars also tend to be in bet^r 
agreement for smaller errors. From our error distribution (Section!^, many sources tend to straddle the limit of Sa reliability. These stars 
also represent the faintest sources, and thus suffer from small numoer statistics when matched to M04. 
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Figu re 1 0. S ame as Figure!^ but residuals with M14. Equa¬ 
tions ( |11[ | and ( |12[ | approximately trace the same levels of reliabil¬ 
ity for m 04 ana ivil4, with reliability increasing for fainter, redder 
stars (r — 2 > 3). 


proper motion errors than we did for M04. As stated 
above, cutting on proper motion errors can yield higher 
reliability. 

Our comparison with LSPM shows us that our cata¬ 
log is reliable for the brightest and faster moving stars. 
The majority of stars with discrepant proper motions 
were found to have more reliable measurements in our 
catalog. Comparisons with M04 and M14 show that our 
catalog tends to be more reliable for brighter sources, but 
reliability remains high for fainter, redder stars. Rather 
than make further cuts on our catalog, we summarize our 
findings and use the results of our comparisons to suggest 
selection criteria to obtain a clean sample in Section 

5. AUGMENTING OUR PROPER MOTIONS WITH 
SDSS-I-USNO-B 


then the nearest matching USNO-B object itself 
matched more than 1 SDSS object. 

2. SIGRA and SIGDEC, the RMS residuals for the 
proper motion fit (in R.A. and Dec., respectively). 

3. nFit, the number of detections used in the fit in¬ 
cluding the SDSS detection (thus, the number of 
plates the object was detected on in USNO-B plus 
one). 

4. O and J, the recalibrated USNO-B O and J mag¬ 
nitudes, respectively, recalibrated to SDSS g. 


We chose to keep the criterion dist22 > 7 since this 
criterion is common among all methods for selecting a 
“clean” sample. We explored each of the criterion out¬ 
lined in the aforementioned studies, seeing how proper 
motion agreement changed with each criteria. 

To select a clean sample, we chose a subset of our color- 
magnitude space where the absolute residuals with M04 
were typically < 10 mas yr“^. The region we selected 
from is shown in Figure In this region, we had 98% 
agreement between M04 and our IF/5'i?-|-SDSS-|-2MASS 
proper motions at 2a (with 381,188 stars). Next, we 
performed a match between our catalog to M04 with 
no requirement on any of the above parameters. This 
match gave us 495,385 stars within the aforementioned 
color-magnitude space. We performed a number of tests, 
using criteria recommended in the various papers listed 
above, and compared the 2cr agreement between ours 
and the M04 proper motion components. T he results of 
our comparisons are summarized in Figure 12 and our 
recommended criteria for selecting high reliability candi¬ 
dates is: 


1. dist 22 > 7 and, 

2. SIGRA < 525 and SIGDEC < 525 and. 


Our combined IF75'E-|-SDSS-|-2MASS proper motions 
are mostly sensitive to faster moving (/Xtot > 20 mas 
yr“^) and redder stars. To make the most complete set 
of proper motions and to ensure the highest completeness 
for low-mass stars within the SDSS photometric sample, 
we require a baseline more sensitive to slower moving 
(tangentially) stars, or disk stars. The ideal remedy is 
M04, which has a higher precision < 4 mas yr”!}, 

but is more sensitive to brighter sources (see Section®. 
To augment our derived proper motions with those from 
M04, we must first choose criteria that select the most 
reliable proper motions from M04. A number of differ¬ 
ent criteria have been proposed for s electing a “clean' 


Dhital et al. 

2010 

Dong et al. 

2011 

West et al. 2011). 

We performed an in-depth exp 
different criteria on selecting a 
The most important parame 
alog for our investigation were 
Server^° ProperMotions Gas 

1. MATCH, the number of c 
matched this object withh 

http://skyserver.sdss 

oration into the effect of 
clean sample. 

;ers within the M04 cat- 
(taken from SDSS Sky- 
Jobs table): 

)bjects in USNO-B that 
1 a 1" radius. If negative, 

org/dr12/en/help/browser/ 


browser.aspx#&&:history=description+ProperMotions+U 


3. (nFit = 6 and MATCH > 1) or 

4. (nFit = 5 and MATCH = 1 and [O < 2 or J < 2]). 

Using these criteria should yield proper motions with 
more than 95% confidence. Another optional cut is to 
remove (redder) objects with r — z > 3, which typically 
have less reliable proper motions, however, these objects 
make up only a small fraction of the objects within M04 
(< 0 . 2 %). 

We chose to supplement our proper motions with those 
from M04, applying all of the above criteria outlined, and 
also requiring a minimum total proper motion greater 
than two times the combined total proper motion error 
ifJ-tot > 2(T^t„J. This gave us 6,620,838 stars, 5,216,855 
of which we did not have a prior measurement for. The 
final number of stars in our catalog for each subsample 
is listed in Table [H 

6 . THE MoVeRS CATALOG 

Our goal was not to construct a “complete” proper mo¬ 
tion catalog, rather, we have attempted to build a cat¬ 
alog of low-mass stars with high-fidelity proper motion 
measurements. We expect our catalog to be more “com¬ 
plete” at the red end of the main-sequence than most 
previous proper motion catalogs due to the use of deeper 
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WISE+SBSS+2MASS 



SDSS+2MASS 



WISE+SBSS 



Figure 11. Same as Figure]^ but for M14 matches. Top: The IF/S'F+SDSS+2MASS stars tend to have more scatter in their residuals 
than the M04 stars, however, the majority are in 3cr agreement with M14. Middle: The SDSS+2MASS stars tend to be in better agreement 
for smaller errors, similar to M04 matches. Bottom: The WISE-\-SDSS stars also tend to be in better agreement for smaller errors, similar 
to M04 matches. The inset plot shows only stars with r — z > 3, which are typically more reliable. 
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NFIT= 4 

NFIT= 5 

NFIT= 6 

525 < (SIGRa & SIGDec) < 1000 

O < 2 or J < 2 

MATCH > 1 

NFIT= 4 

94.9% 

(19279) 






NFIT= 5 


96.1% 

(88247) 





NFIT= 6 



98.1% 

(381188) 




525 < (SIGRA & SIGDEC) < 1000 

35.5% 

(34) 

39.3% 

(89) 

75.5% 

(204) 

75.5% 

(204) 



O < 2 or J < 2 

93.1% 

(13340) 

96.0% 

(66019) 

96.1% 

(538) 

50.0% 

(2) 

97.3% 

(518) 


MATCH > 1 

88.3% 

(571) 

88.7% 

(477) 

97.9% 

(662) 

100.0% 

(2) 

100.0% 

(5) 

98.6% 

(639) 


Figure 12. Table showing the percentage of stars that had a 2c7 agreement in both proper motion components {a and (5) between our 
IT75'F+SDSS+2MASS proper motions and those from M04, as a function of M04 parameters. The number of stars that met each criteria 
are listed below the percentage in parenthesis. Colors represent our limits for selecting: 1) “good” criteria (green); 2) “bad” criteria (red); 
3) baseline, or starting, criteria (blue); and 4) criteria that did not include enough stars to be meaningful (light gray). Default criteria, 
unless changed, were: 1) nFit = 6; 2) SIGRA & SIGDEC < 525; and 3) MATCH = 1. 
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Figure 13. r — z distributions for each subsample in our cat alog. 
Approximate spectral types taken from Hawley et al.| |2002| and 
|Bochanski et al.|(|2007||. M04 and our SL!bb+2MAbb distributions 
typically include more bluer stars, while proper motions measured 
with WISE typically include more redder stars. 


surveys, specifically in the IR. To illustrate this point, 
we show the r — z color distribution for each of our sub¬ 
samples (e.g., W75'i?-I-SDSS-I-2MASS or SDSS-I-USNO- 
B) in Figure 13 Our W75'i?-|-SDSS-|-2MASS baseline 


contains a large number of star with r — z > 2.5, and our 
WISE+SDSS subsample also contains more stars with 
r — z > 2.5 than our M04 stars, which peak at a spectral 
type of ^M4. Our newly computed proper motions will 
identify many new very-low-mass stars with r — z > 2.5 
(spectral types later than ^M6). The schema for our 
catalog can be found in Appendix 0 Proper motions 
and errors for ea ch o f the subsamples in our catalog are 
shown in Figure 14 The extremely wide wings for our 
WTS'A-I-SDSS anaSDSS-|-2MASS matches are due pri¬ 
marily to stars with short time baselines (< 3 years), 
which typically overestimate proper motions. 

A small number of our stars have high proper mo¬ 
tions {utot > 1000 mas yr“^; ~0.5%). Some of these 
stars are due to our large search radius (6"), which may 
pull in neighboring objects. Our color selection criteria 
should remove a number of these, however, in crowded 
fields (e.g., within the Galactic plane), these objects are 
more prominent. Rather than remove fast moving ob¬ 
jects, since many of them will be true detections, we 
recommend a conservative cut to remove potential out¬ 
liers is to eliminate stars near the Galactic plane (e.g.. 


| 6 | < 20 °). 

Another potential cut is to remove objects that have 
appeared to move a distance close to our search radius 
over the time baseline (e.g., objects that have moved ^6" 
within a 1 year period). We plan to determine the valid¬ 
ity of these high proper motions sources in a future study. 
Also, selecting stars that either satisfy equation (111 
or (12) can yield higher reliability proper motions, but 
is biased towards selecting stars with all three epochs. 
Lastly, selecting stars with smaller proper motion errors 
(cr^tot ^ 60 mas yr“^), specifically for the 1F75A-I-SDSS 
and SDSS-(-2MASS subsamples can increase reliability. 
The M04 


7. DISCUSSION 

This paper presents an improved proper motion cat¬ 
alog for low-mass stars, specihcally, completing the red 
end of the low-mass sequence. Based off our color se¬ 
lection, this catalog spans from late K dwarfs to early 
L dwarfs. In particular, this catalog will allow unprece¬ 
dented studies of subdwarfs and potential halo stars. To 
show the relative populations of stars within o ur cata¬ 
log, w e include reduced proper motions (RPMs; Luyten 
19221), given as 


Hr = r + 5 + 51og^ = Mr — 3.25 -|- 51ogTT) (13) 

where /r is the total proper motion in arcsec yr“^, and 
vt is the heliocentric tangential velocity in km s“^ given 
by Vt = 4.74 X {fj, ■ d), where d is the distance in par¬ 
secs. RPM diagrams for each of our four subsamples 
are shown in Figure]^ To segregate different kinematic 
groups, we have drawn a line at vt = 180 km s“^, the 
approximat e limit which dif ferentiates disk stars from 
subdwarfs ( Sesar et al.| 2008). Although subdwarfs may 
scatter above the 180 km s“^ line, disk d warfs are not 
typic ally found below the 180 km s“^ line (Dhital et al. 

1oio|. 


The SDSS-bUSNO-B stars are primarily earlier type 
disk stars. Our 1F/S'i?-|-SDSS-|-2MASS stars are also pri¬ 
marily disk dwarfs, but peak at redder colors, and are 
more complete (and reliable) at the reddest end of the 
main-sequence. Our SDSS-I-2MASS and WTS'A-I-SDSS 
baselines appear to probe bluer and redder subdwarf 
populations, respectively. 

We have identified a number of science questions which 
can be investigated with this catalog: 

1. Identify SDSS spectroscopic high-proper motion 
low-mass stars to find hypervelocity candidates 
(Favia et al., submitted.). 

2. Identify extremely low-mass, common proper mo¬ 
tion binaries. 


3. Investigate Galactic kinematics for the lowest-mass 
main-sequence members. 

4. Identify low-mass field stars with infrared excesses. 

5. Confirm pervious catalogs of wide-binaries (e.g., 
Dhital et aL]|2015 ). 

With current and future efforts multi-epoch surveys, such 
as Gaia and LSST, this catalog will also prove invaluable 
as a calibrator for the reddest stellar populations within 
these surveys. 
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Figure 14. Proper motion and proper motion error distributions for our catalog. The vast majority of our catalog is comprised of stars 
with fitot < 1000 mas yr”'^ (99.5% of the entire catalog). 


7.1. Common Proper Motion Binaries: Investigating 
the SloWPoKES-II Catalog 

To demonstrate the utility of MoVeRS, we matched our 
catalog to the second release of the Sloan Low-mass Wide 
Pa irs of Kinematicall y Equivalent Stars (SLoWPoKES- 
11; Dhital et al. 2015), a sample of wide binaries identi- 
fied without prope r motions. Th is catalog is an ex tension 
of the study from |Dhital et al.| (SLoWPoKES-I; |2010|), 
where it was shown that wide-binaries with separations 
less than ^20" could be identified based off similar dis¬ 
tances, but without the need for proper motions. There 
were 260 matches between the two catalogs. Figure 16 
shows the distributions for the proper motion compo¬ 
nents between the prima ry and secondary, t he distribu¬ 
tion of equation (6) from Dhital et al. (2010), a measure 
of the weighted difference in both proper motion compo¬ 
nents given as 


Xr. = 


f A/To 


\ ^Afio 


+ 




(14) 


The fraction of matches that met the criteria of reli¬ 


able proper motions binaries < 2; Dhital et al. 


2010) was 38%. If we restrict this analysis to pairs 


that have proper motions with all three epochs (i.e. 
1T/5E-I-SDSS-I-2MASS; 127 pairs), this fraction increase 
to 47%. 

The original SLoWPoKES catalog did not have a for¬ 
mal constraint on average distance due to the fact that 
the range of spectral types (mid-K to mid-M) made 
the catalog sensitivity a function of color and distance. 
It was found that reliability was higher for stars with 
d < 1200 pc, due to smaller photometric distance un¬ 
certainties. Since SLoWPoKES-11 has no constraint on 
average distance, we chose to explor e h ow pair fidelity 
corresponded with distance. Figure [17] shows the dis¬ 
tribution of both high-fidelity pairs (^ < 2) and low- 
fidelity pairs (xti > 2). High-fidelity pairs are primarily 
found within 1200 pc, with fidelity becoming noticeably 
worse at larger distances. 

We chose to investigate how angular separation be¬ 
tween the primary and secondary component affected the 
likelihood of pairs having common proper motions. From 
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SDSS+USNO-B VE/SE+SDSS+2MASS SDSS+2MASS WISE+SBSS 
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Figure 15. Reduced proper motion diagrams for the four sub-samples included in our c atalog, each bin is (0.1 mags)^. The dashed line 
represent a tangential velocity of 180 km s“^, which separates disk stars from halo stars (|Sesar et al.|2008ll. Our catalog provides one of 
the largest samples of high-confidence subdwarfs (stars below the dashed line). The M04 source at r — z Ri 2, Hr ~ 30 appears to have a 
spurious proper motion (as observed in archived images). 


the results of Dhital et al. 


(2010) and Dhital et al. (2015), prove the fidelity of the sample. 


we expect closer pairs to have a higher probability ot be- 
ing true binaries. However, we only have the ability to 
probe pairs with separations > 8", except for a small sub¬ 
set of objects matched in SDSS+2MASS. To investigate 
reliability as a function of pair separation, we calculated 
the fraction of reliable pairs in a moving l"bin and com¬ 
puted binomial uncertainties, using only pairs with an 
average distance < 1200 pc since reliability was shown 
to decrease significantly past this point. Our results are 
shown in Figure [TS} as expected, proper motion reliabil¬ 
ity decreases with angular separation; for sources with 
angular separations > 15", we find reliability drops to 
zero. Due to the small sample size of our matched pairs, 
we do not attempt to make any claims on the overall reli¬ 
ability of SLoWPoKES-II; this exercise was simply used 
to show the applicability of our catalog. From our ini¬ 
tial results, coupled with the fact that the distribution of 
angular separations in the SLoWPoKES-II catalog peaks 
at separations < 8", we expect many of the pairs to be 
bona-fide binaries. Choosing pairs with smaller average 
distances and smaller angular separations can help im- 


8 . SUMMARY 

We have created a catalog containing 8,735,004 proper 
motion verified photometric low-mass stars. Proper mo¬ 
tions were computed using the WISE, SDSS, and 2MASS 
surveys (3,518,150 stars), and augmented with proper 
motions from SD SS+USNO-B (5,216,854 stars; iMunn 


et al.|2004 2008). All stars were required to have a total 


proper motion greater than twice the uncertainty in their 
measurement, thus ensuring high-fidelity main-sequence 
stars. The estimated precision of our catalog is ^10 
mas yr“^ for our 1F/S'E+SDSS-|-2MASS sources, and 
^40 mas yr“^ for our SDSS-I-2MASS and WTS'E+SDSS 
sources, primarily due to shorter time baselines. 

Comparison against the high proper motions stars from 
LSPM suggests good agreement for high proper mo¬ 
tion stars in our catalog in all three subsamples (e.g., 
WTS'E-I-SDSS or SDSS+2MASS). For our subsamples, 
agreement with LSPM at the 2cr level was 98%, 97% 
and 96% for lF/5E-bSDSS+2MASS, SDSS-b2MASS, and 
WISE-\-SDSS, respectively (these all increase to 99% at 
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Figure 16. Distributions for proper motion residuals between the primary and secondary components for the stars matched between our 
MoVeRS catalog and SLoWPoKES-IL The solid line denotes all pairs matched, and the dotted line denotes only stars matched with all 
three epochs (i.e. WISE-\-SDSS-\-2MASS). Subplots (a) and (b) show the proper motion differences between the primary and secondary in 
a and 5, respectivel y. B oth distributions are peaked at zero (dashed lines), showing good agreement between the proper motions. Subplot 
(c) shows equation (|14[ ), the quadrature sum of the weighted difference of both proper motion components. The red dotted line denotes a 
value of 2, the cuton value for binaries to be include in the original SLoWPoKES catalog. 
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Figure 17. Average distance distributions for pairs from the 
SLoWPoKES-II catalog. Dark gray lines correspond to high- 
fidelity pairs (X(i < 2; equation ITS, and light grey to low-fidelity 
{Xfj. > 2) pairs. Dotted lines correspond to all pairs found within 
MoVeRS, while solid lines correspond to only pairs where both 
components had three epochs (i.e. IT75'i?-|-SDSS+2MASS). Pairs 
with an average distance > 1200 pc (dashed line) tend to have less 
similar proper motions (and hence lower fidelity) than closer pairs. 
The majority of pairs within SLoWPoKES-II have angular separa¬ 
tions < 8”, so we expect many of the pairs to have high-reliability. 


the 3(7 level). We further compared our proper motions 
to SDSS-hUSNO-B measur ements and the deep er proper 
motion catalog released by Munn et al. (2014|). In both 
cases, our proper motion precision is strongly correlated 
with apparent magnitude and color, diminishing for bluer 
and fainter sources. 

The utility of this catalog will be in the vast number of 
motion verified low-mass stars it contains, and its high 
reliability, specifically for the reddest and lowest-mass 
members of the catalog. We expect the red end of this 
catalog to surpass the limits of Gaia. Our catalog is 
available through SDSS CasJobs and VizieR. 
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Table 5 

i — z Color Selection Criteria 


r — i 

i — Z 

^ of Stars 

Bin Size 

0.200 

0.367 ± 0.044 

309784 

0.01 

0.205 

0.367 ± 0.044 

309784 

0.01 

0.215 

0.379 ± 0.048 

318938 

0.01 

2.010 

2.685 ± 0.210 

8 

0.2 

2.020 

2.685 ± 0.210 

8 

0.2 

2.030 

2.685 ± 0.210 

8 

0.2 


Note. — Table is published in its en¬ 
tirety in the electronic edition of AJ, a portion 
is shown here for guidance regarding its form 
and content. 
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Collaboration et al. 2013). Figures in this work were 
created usi ng the Fytho ii based graphics environment 
Matplotlib (Hunter||2007). 


Figure 18. Reliability fraction of pairs (with reliable pairs defined as pairs with x^^ < 2) as a function of angular separation between 
components. The bin size is shown in the top left corner. Only pairs with average distances < 1200 pc were used. Reliability fractions were 
calculated for all matched pairs (light gray square), and pairs where both components had inSi?+SDSS+2MASS measurements (dark 
gray circles). We did not find any reliable pairs with angular separations > 15”. Reliability increases for smaller angular separations, 
unfortunately we cannot probe the majority of the SLoWPoKES-II pai rs, which have anefular separatio ns < 8” 

Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro¬ 
cessing and Analysis Center/California Institute of Tech¬ 
nology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 

This publication also makes use of data products from 
the Wide-field Infrared Survey Explorer, which is a joint 
project of the University of California, Los Angeles, 
and the Jet Propulsion Laboratory/California Institute 
of Technology, funded by the National Aeronautics and 
Space Administration. 

This research made use of Astropy, a co mmunity- 


The authors are also pleased to acknowledge that much 
of the computational work reported on in this paper was 
performed on the Shared Computing Cluster which is ad¬ 
ministered by Boston University’s Research Computing 
Services (www.bu.edu/tech/support/research/). In 
particular, C.A.'i'. would like to thank Paul Ualba and 
Yann Tambouret for helpful discussions regarding opti¬ 
mization of code and parallelization. 


developed core Python package for Astronomy (Astropy 


APPENDIX 

A.l COLOR SELECTION POLYGONS 

Tables EE and contain the color selection criteria we used to trace the stellar locus and select our initial sample. 


A.2 QUERYING THE CATALOG 

Our Motion Verified Red Stars (MoVeRS) catalog is available through SDSS CasJobs^^ and VizieR^^. To access our 
catalog through CasJobs, please refer to the documentation for accessing public tables. The following is an example 
SQL query for accessing our table within the DRIO context to return SDSS positions, rizJHKWl photometry, proper 
motions, and proper motion errors for stars with: 1) total proper motions less then 500 mas yr“^; 2) total proper 
motion errors less than 20 mas yr“^; 3) r — z > 2.5; and 4) 18 < r < 21: 


http://skyserver.sdss.org/casiobs/ 

http://vizier.u-strasbg.fr/viz-bin/VizieR 
























22 


Theissen, West, & Dhital 


Table 6 

z — J Color Selection Criteria 


r — 2 

z — J 

of Stars 

Bin Size 

0.500 

1.030 ±0.056 

1384 

0.01 

0.505 

1.030 ±0.056 

1384 

0.01 

0.515 

1.033 ±0.054 

3861 

0.01 

4.900 

2.076 ± 0.049 

3 

0.2 

4.910 

2.044 ± 0.024 

2 

0.2 

4.920 

2.044 ± 0.024 

2 

0.2 


Note. — Tablej^is published in its entirety in 
the electronic edition of AJ, a portion is shown 
here for guidance regarding its form and content. 


Table 7 

z — W1 Color Selection Criteria 


r — 2 

z-Wl 

of Stars 

Bin Size 

0.500 

1.712 ± 0.096 

4889 

0.01 

0.505 

1.712 ± 0.096 

4889 

0.01 

0.515 

1.724 ±0.096 

13629 

0.01 

4.900 

3.429 ± 0.082 

3 

0.2 

4.910 

3.373 ±0.027 

2 

0.2 

4.920 

3.373 ±0.027 

2 

0.2 


Note. — Tablej^is published in its entirety in 
the electronic edition of AJ, a portion is shown 
here for guidance regarding its form and content. 


SELECT 

p.sdss_ra, p.sdss_dec, p.rmag, p.imag, p.zmag, p.jmag, p.hmag, p.kmag, 
p.wlmpro, p.pmra, p.pmdec, p.pmra_toterr, p.pindec_toterr 
FROM public.LowMassPM.MoVeRS p 
WHERE 

p.pmra * p.pmra + p.pmdec * p.pmdec < 500 * 500 AND 

p.pmra_toterr * p.pmra_toterr + p.pmdec_toterr * p.pmdec_toterr < 20 * 20 AND 
p.rmag - p.zmag > 2.5 AND 
p.rmag BETWEEN 18 and 21 


A.3 MOVERS CATALOG SCHEMA 


Table 8 
Catalog Schema 


Field Name Format Units Description 


SDSS_OBJID 

int64 


SDSS_RA 

float 64 

degrees 

SDSS_DEC 

float64 

degrees 

SDSS_RAERR 

float 3 2 

degrees 

SDSS_DECERR 

float 3 2 

degrees 

SDSS_MJD 

float 3 2 

days 

UMAG 

float 3 2 

mags 

GMAG 

float 3 2 

mags 

RMAG 

float 3 2 

mags 

IMAG 

float 3 2 

mags 

ZMAG 

float 3 2 

mags 

UMAG_ERR 

float 3 2 

mags 

GMAG_ERR 

float 3 2 

mags 

RMAG_ERR 

float 3 2 

mags 

IMAG_ERR 

float 3 2 

mags 

ZMAG_ERR 

float 3 2 

mags 

2MASS_RA 

float 3 2 

degrees 

2MASS_DEG 

float 3 2 

degrees 

2MASS_RAERR 

float 3 2 

degrees 

2MASS_DEGERR 

float 3 2 

degrees 

2MASS_MJD 

float 3 2 

days 

2MASS_PH_QUAL 

3 character string 


2MASS_RD_FLG 

3 character string 


2MASS BL±'LG 

3 character string 



SDSS DR8+ Object ID 
SDSS R.A. 

SDSS Decl. 

SDSS R.A. error (in proper units, i.e. Aacos J) 
SDSS Decl. error 

SDSS r-band modified Julian date 
SDSS u-band PSF magnitude 
SDSS 5'-band PSF magnitude 
SDSS r-band PSF magnitude 
SDSS i-band PSF magnitude 
SDSS 2-band PSF magnitude 
SDSS u-band PSF magnitude error 
SDSS p-band PSF magnitude error 
SDSS r-band PSF magnitude error 
SDSS 2-band PSF magnitude error 
SDSS 2-band PSF magnitude error 
2MASS R.A. 

2MASS Decl. 

2MASS R.A. error (in proper units, i.e. AqcosJ) 

2MASS Decl. error 

2MASS modified Julian date 

2MASS photometric quality fiag 

2MASS read flag 

2MASS blend flag 
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Table 8 — Continued 


Field Name 

Format 

Units 

Description 

2MASS_CC_FLG 

3 character string 


2MASS contamination and confusion flag 

2MASS_GAL_CONTAM 

int32 


2MASS extended source “contamination” flag 

JMAG 

float 3 2 

mags 

2MASS J-band PSF magnitude 

JMAG_ERR 

float 3 2 

mags 

2MASS J-band PSF corrected magnitude uncertainty 

JMAG_ERRTOT 

float 3 2 

mags 

2MASS J-band PSF total magnitude uncertainty 

JSNR 

float 3 2 

mags 

2MASS J-band SNR 

HMAG 

float 3 2 

mags 

2MASS iJ-band PSF magnitude 

HMAG_ERR 

float 3 2 

mags 

2MASS iJ-band PSF corrected magnitude uncertainty 

HMAG_ERRTOT 

float 3 2 

mags 

2MASS JJ-band PSF total magnitude uncertainty 

HSNR 

float 3 2 

mags 

2MASS H-band SNR 

KMAG 

float 3 2 

mags 

2MASS A's-band PSF magnitude 

KMAG_ERR 

float 3 2 

mags 

2MASS ATs-band PSF corrected magnitude uncertainty 

KMAG-ERRTOT 

float 3 2 

mags 

2MASS ATs-band PSF total magnitude uncertainty 

KSNR 

float 3 2 

mags 

2MASS RTs-band SNR 

J_PSECHI 

float 3 2 


2MASS J-band reduced goodness-of-fit for the PSF 

H_PSFCHI 

float 3 2 


2MASS iJ-band reduced x^ goodness-of-fit for the PSF 

K_PSFCHI 

float 3 2 


2MASS ifs-band reduced x^ goodness-of-fit for the PSF 

WISE_RA 

float 3 2 

degrees 

WISE R.A. 

WISE_DEC 

float 3 2 

degrees 

WISE Decl. 

WISE_RAERR 

float 3 2 

degrees 

WISE R.A. error (in proper units, i.e. Aq;cos(5) 

WISE_DECERR 

float 3 2 

degrees 

WISE Decl. error 

WISE_GC_FLG 

4 character string 


WISE contamination and confusion flag 

WISE_EXT_FLG 

int32 


WISE extended source flag 

WISE_VAR_FLG 

4 character string 


WISE variability flag 

WISE_PH_QUAL 

4 character string 


WISE photometric quality flag 

WISE_W1MJDMEAN 

float 3 2 

days 

WISE VPl-band average modified Julian date 

WIMJDSIG 

float 3 2 

days 

WISE MJD uncertainty'^ 

WIMPRO 

float 64 

mags 

WISE VFl-band PSF magnitude 

WISIGMPRO 

float64 

mags 

WISE VFl-band PSF magnitude uncertainty 

WISNR 

float 64 


WISE VFl-band SNR 

W1RCHI2 

float 3 2 


WISE reduced x^ goodness-of-fit for the PSF 

NEAREST_NEIGHBOR 

float 3 2 

arcsec 

Distance to nearest SDSS primary object 

NEAREST_RMAG 

float 3 2 

mags 

SDSS r-band PSF magnitude of nearest neighbor 

NEAREST JMAG 

float 3 2 

mags 

SDSS 2 -band PSF magnitude of nearest neighbor 

NEAREST_ZMAG 

float 3 2 

mags 

SDSS 2 :-band PSF magnitude of nearest neighbor 

NEIGHBORS 

int32 


Number of SDSS primary objects within 15^^ 

RRl 

int32 


Flag if there is an object within 8^^ with Tgource — ^"neighbor > — 1 

RR2 

int32 


Flag if there is an obiect within 8^^ with rsmime — > —2 

RR25 

int32 


Flag if there is an obiect within 8^^ with Tsmirre — > —2.5 

RR3 

int32 


Flag if there is an object within 8'' with Tgource — ^"neighbor ^ 

RR4 

int32 


Flag if there is an object within 8^^ with rgource — "Tneighbor ^ 

RR5 

int32 


Flag if there is an object within 8^^ with rgource — rneighbor ^ 

PMRA 

float 3 2 

mas yr~^ 

Proper motion in R.A. (in proper units, i.e. /ia cos J) 

PMDEC 

float 3 2 

mas yr~^ 

Proper motion in Decl. 

PMRA_M04 

float 3 2 

mas yr~^ 

M04 Proper motion in R.A. (in proper units, i.e. fia cosJ) 

PMDEC_M04 

float 3 2 

mas yr~^ 

M04 Proper motion in Decl. 

PMRAJNTERR 

float 3 2 

mas yr~^ 

Intrinsic error in proper motion in R.A. 

PMDEC JNTERR 

float 3 2 

mas yr~^ 

Intrinsic error in proper motion in Decl. 

PMRA_MEASERR 

float 3 2 

mas yr~^ 

Measurement error in proper motion in R.A. 

PMDEC-MEASERR 

float 3 2 

mas yr~^ 

Measurement error in proper motion in Decl. 

PMRA_FITERR 

float 3 2 

mas yr~^ 

Fit error in proper motion in R.A. 

PMDEC_FITERR 

float 3 2 

mas yr~^ 

Fit error in proper motion in Decl. 

PMRA.TOTERR 

float 3 2 

mas yr~^ 

Combined error in proper motion in R.A. 

PMDEC.TOTERR 

float 3 2 

mas yr~^ 

Combined error in proper motion in Decl. 

PMRAERR_M04 

float 3 2 

mas yr~^ 

M04 error in proper motion in R.A. 

PMDECERR_M04 

float 3 2 

mas yr~^ 

M04 error in proper motion in Decl. 

BASELINE 

float 3 2 

years 

Time baseline used to compute our proper motions 

DBIT 

3 character string 


Detection bit identifving survevs used in computing proper motions^ 

RECOMP 

int32 


Flag indicating proper motions were recomputed (see Section|3.4.l|l 

USE 

int32 


Flag indicating which PM measurement to use^ 

MATCH_M04 

int32 


Number of SDSS objects within a 1'^ radius matching the USNO-B object 

SIGRA_M04 

float 3 2 

mas 

M04 RMS residual for the proper motion fit in R.A. 

SIGDEC_M04 

float 3 2 

mas 

M04 RMS residual for the proper motion fit in Decl. 

NEIT_M04 

int32 


Number of detections used in the M04 fit 

O_M04 

float 3 2 

mags 

Recalibrated USNO-B O magnitude, recalibrated to SDSS g 

J_M04 

float 3 2 

mags 

Recalibrated USNO-B J magnitude, recalibrated to SDSS g 

WS_DIST 

float 3 2 

arcsec 

Total distance between WISE position and SDSS position^ 

S2_DIST 

float 3 2 

arcsec 

Total distance between SDSS position and 2MASS position ^ 

W2_DIST 

float 3 2 

arcsec 

Total distance between WISE position and 2MASS position^ 
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Table 8 — Continued 


Field Name Format Units Description 


Defined as .5x (WIMJDMAX-WIMJDMIN). 

‘111’: WISE, SDSS, and 2MASS were used; ‘110’: SDSS and 2MASS were used; ‘Oil’: WISE and SDSS were used; ‘000’: SDSS+USNO- 
B measurement is available. 

^ ‘1’: proper motions were measured here; ‘2’: proper motions are from M04; or ‘3’: both proper motions are available. 

^ Total distance = ^/(AQp~cos'^i~cos'^2~+T^^T^- 
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